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I. INTRODUCTION 



Hard exclusive reactions and generalized parton distributions (GPDs) have been in the 
focus of hadronic physics for the last decade GPDs interpolate between 

elastic form factors and structure functions and contain detailed information on distributions 
and correlations of partons (quarks and gluons) in hadronic targets (pions, nucleons, and 
nuclei). In particular, GPDs describe the distribution of partons both in the longitudinal 
momentum direction and in the impact parameter (transverse) plane [3] and also allow to 
access the total angular momentum of the target carried by the partons jsj]. 

The QCD factorization theorems for hard exclusive processes 9, [lo| state that GPDs are 
universal distributions that enter the perturbative QCD description of various hard exclusive 
processes such as Deeply Virtual Compton scattering (DVCS), 7* T — ► 7T (T denotes 
any hadronic target), exclusive production of mesons, 7* T — > MT [where M denotes a 
(pesudo) scalar or a vector meson], and many other processes, including generalizations of 
these two reactions. 

Although the factorization theorems make it theoretically possible to extract GPDs from 
the data, this is a difficult task in practice since GPDs are functions of four variables and 
the GPDs enter experimentally measured observables in the form of convolution with hard 
coefficient functions. Therefore, there is a clear need for modeling GPDs, both to interpet 
the results of the completed experiments in terms of the microscopic structure of the hadron 
target and also to plan future experiments. 

In this work, we study quark and gluon GPDs of heavy nuclei and DVCS on nuclear 
targets at small values of Bjorken xb (large energies) . In particular, we generalize the 



theory of leading twist nuclear shadowing 



11 



12 



13] to the case of GPDs and compute 
next-to-leading order quark and gluon GPDs of nuclei for 10~ 5 < ig < 0.2 and at a fixed 
virtuality Q 2 . Using the obtained nuclear GPDs, we compute the DVCS amplitude, the 
DVCS cross section, and the DVCS beam-spin asymmetry for the heavy nuclear target of 



208 



Pb. Our results can be summarized as follows: 

(i) Leading twist nuclear shadowing suppresses very significantly the DVCS amplitude 
and the DVCS cross section at small values of Bjorken xb- 

(ii) In the £ — > limit, nuclear GPDs reduce to impact-parameter-dependent nuclear 
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parton distribution functions (PDFs). Therefore, nuclear GPDs allow one to access 
the spatial image of nuclear shadowing. The shadowing correction to nuclear GPDs 
introduces nontrivial correlations between the light-cone fraction x and the impact 
parameter b. 

(iii) DVCS interferes with the purely electromagnetic Bethe-Heitler (BH) process. At small 
values of the momentum transfer t, which dominate coherent nuclear DVCS (without 
nuclear break-up), and also for the t-integrated cross sections, the BH cross section is 
much larger than the DVCS one. This makes it rather challenging to extract a small 
DVCS signal on the background of the dominant BH contribution to the eA — > ejA 
cross section. However, owing to the rapid t-dependence, the DVCS cross section 
becomes (much) larger than the BH cross section near the minima of the nuclear 
form factor. This suggests that the measurements of nuclear DVCS at the values of t 
close to the minima of the nuclear form factor will not only be very sensitive to the 
magnitude of nuclear shadowing (owing to the suppression of the nonshadowed Born 
contribution), but will also have a sufficiently small Bethe-Heitler contribution. 

(iv) Another possible way to access nuclear GPDs in the small xb region is through the 
measurement of the DVCS beam-spin asymmetry, Ajjj. Nuclear shadowing causes 
dramatic oscillations of the asymmetry at the fixed (f) = 90° as a function of the 
momentum transfer t. The position of the points where Alu changes sign is directly 
related to the magnitude of nuclear shadowing. 

The rest of the paper is structured as follows. In Sec. [Til we derive the expression for 
nuclear shadowing for nuclear GPDs. In Sec. IHH we analyze the £a —> limit of the resulting 
nuclear GPDs, point out the equivalence of the nuclear GPDs in this limit to the impact- 
parameter-dependent nuclear PDFs, and discuss the spacial image of nuclear shadowing. 
Predictions for DVCS observables (the DVCS amplitude and cross section and the beam- 
spin DVCS asymmetry) are presented and discussed in Sec. IIVI Finally, we summarize and 
draw conclusions in Sec. [Vj 
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II. LEADING TWIST NUCLEAR SHADOWING AND NUCLEAR GPDS 



The nuclear structure function F 2 a(x b ,Q 2 ) measured in inclusive deep inelastic scat- 
tering (DIS) with nuclear targets differs from the sum of free nucleon structure functions 
F 2N {x B ,Q 2 ) over the entire range of values of Bjorken x B , lly, [l7|. In particular, 



for small values of x B , 10~ 5 < x B < 0.05 - 0.1, F 2A (x B , Q 2 )/[AF 2N (x B , Q 2 )} < 1, which is 
called nuclear shadowing. 

As we learned from DIS with fixed nuclear targets, the effect of nuclear shadowing is 
quite large for small x B . The kinematics of the future high-energy collider Q] will cover 



the small-a:^ region, where the effect of nuclear shado wing wil 



11 



12 



play a major role. 
131 ] is an approach to nuclear 



The leading twist (LT) theory of nuclear shadowing 
shadowing, in which nuclear shadowing in DIS with nuclei is explained in terms of hard 
diffraction in lepton- nucleon DIS. In particular, by using the QCD factorization theorems 
for inclusive and hard diffractive DIS and generali zing the result for nuclear shadowing in 
pion-deuteron scattering obtained by V.N. Gribov 201] . the leading twist theory of nuclear 
shadowing makes predictions for the shadowing correction to nuclear PDFs, Sfj/A(x B , Q 2 ) = 
fj/A(x B ,Q 2 ) — Afj/ N (x B ,Q 2 ), in terms of the free nucleon (proton) diffractive PDFs 
for small values of x B , 10~ 5 < x B < 0.2. One should note that the generalization of Gribov's 
result to DIS and to nuclei heavier than deuterium makes an explicit assumption that the 
diffractive state produced in the interaction with the first nucleon of the t arg et elastically 
rescatters off the rest of the nucleons (quasi-eikonal approximation) H, I12L In the 



limit of low nuclear density, when the interaction with only two nucleons of the target is 
important, the relation between 5fj/A(x B ,Q 2 ) and is model independent. Since f^m 
is a leading twist quantity, so is 5fj/A(x B , Q 2 ), which explains the name leading twist theory 
of nuclear shadowing. 

In this work, we generalize the theory of leading twist nuclear shadowing of usual nuclear 
PDFs 0, Q, Q to the off-forward kinematics, DVCS on nuclear targets, and nuclear 
GPDs. The DVCS amplitude on any hadronic target is defined as a matrix element of the 
T-product of two electromagnetic currents (see, e.g., Ref. [3j), 



H% = -1 \ d*xe-^ x (P A \T{J»(x)r(0)}\P A ) 



(1) 



where q (—q 2 = Q 2 ) is the momentum of the virtual photon and Pa and P' A are the momenta 
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of the initial and final nucleus, respectively. DVCS on a nuclear target is presented in Fig. CD 

For the analysis of the matrix element in Eq. (jTJ) , it is convenient to introduce two light-like 
vectors p = l/\/2(l, 0, 0, 1) and n = 1/ V / 2(l, 0, 0, —1) and to work in the so-called symmetric 
frame, where q and the average momentum of the initial and final nucleus, Pa = (Pa+Pa)/2, 
are large and have no transverse component (with respect to the light-like directions defined 
by p and n). Then, the involved momenta can be parameterized as jj] 

jCf 2 X 

Pa = (l + U)PiP+^f(l-U)n-^, 

M 2 A 
P'a = (l-U)PlP+^(l + U)n + ^, 

M 2 

A = P> A -P A = -2UPlP + U 1 £n+ A ± , 
Q = -2UPiP+T~+n, (2) 

where Pj[ = Pa ■ n; M\ = M\ — t/4, with Ma the mass of the nucleus and t = A 2 
the momentum transfer squared; Q 2 is the photon virtuality; A^ is the component of A 
orthogonal to the vectors p and n. As follows from the decomposion of Eq. ([21, 

where xa is the Bjorken variable with respect to the nuclear target, 

XA= 2PA— q = A XB - (4) 
The Bjorken variable xb is defined in the usual way with respect to a free nucleon. 

In this work, we shall consider spinless nuclei since we are not concerned with spin effects 
in nuclear shadowing. To the leading twist accuracy and to the leading order in the strong 
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FIG. 2: Feynman graphs corresponding to the DVCS amplitude on a nuclear target, H A , showing 
the impulse (Born) approximation (a) and the shadowing correction arising from the interaction 
with two nucleons (b) and three nucleons of the target (c), respectively. 



coupling constant, H^ A V of a spinless nucleus is expressed in terms of a single generalized 
parton distribution, Ha-, convoluted with the hard scattering coefficient function C + (x, £4.) 
(see, e.g., Ref. jj), 

= J dx C+(x, U)Ha(x, U, t, Q 2 ) = -ft H A (U, t, Q 2 ) , (5) 

where = g^ u —p^n u —p u n fl ; C + (x, £4) = l/(x — £a + ze) + l/( x + £a — ze). The function 
H.A is also called the Compton form factor (CFF). 

At sufficiently high energies (small Bjorken xb), the virtual photon interacts with many 
(all) nucleons of the target and the DVCS amplitude on a nuclear target, Hjf, receives 
contributions from the graphs presented in Fig. [2j Figures Hl^a), [2(b), and[2|c) correspond 
to the interaction with one, two, and three nucleons, respectively. Graphs that correspond 
to the interaction with four and more nucleons of the target are not shown, but they are 
implied. Therefore, H^ 1 can be written as the following sum: 

Ha = H ( A a)flu + H A b)flu + H A c)flu + . . . , (6) 

where the terms in the right-hand side correspond to the graphs shown in Figs. |2(a), [2(b), 
and [2](c), respectively. 
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A. Impulse approximation 



Let us start with the calculation of the graph shown in Fig. 12(a). For the case of a 



deuterium target, the derivation was done in Ref. 2l|. Therefore, in this subsection, we 
shall follow Ref. 2l| making straightforward generalizations to heavier nuclei and high- 
energy kinematics. 

The calculation of the graph in Fig. 2(a) can be carried out straightforwardly using 
the light-cone (LC) formalism. In this formalism, each state is characterized by its plus- 
momentum, p + = p ■ n = (p° + p 3 )/\^2, the transverse momentum, p±, and the helicity, A. 
The minimal Fock component of the nuclear state \Pa) is expressed in terms of the nuclear 
LC wave function 4>a and the product of nucleon states as 

1 da.i d 2 k\ / [ 1 / 



Aj J i=l V 1 \j=l / \j=l 

x k± u Xi,a 2 , k± 2 , A 2 , . . . )\(XiP][, k ±i + a>iP± A} A») , (7) 



where a.i = pj ' jP\ is the fraction of the nucleus plus-momentum carried by nucleon i. Since 
we are not concerned with the correlations of nucleons in the target nucleus, we take 4>a as 
a product of the light-cone wave functions of individual nucleons, 4>Ni 

A 

4>a{ol\, k±i, Ai, a 2 , k ±2 , A 2 , . . . ) = Jj 0jv(aii, k ±i , A») . (8) 

i=i 

Substituting Eq. (JTj) for the initial and final nuclear states in the nuclear DVCS ampli- 
tude [Eq. fl5J)], we obtain 

x (p' N \T{r(x)r(o)}\ PN } } (9) 

where YIn denotes the sum over active (interacting) nucleons. In Eq. and in the rest 
of the paper, we neglect the off-shellness of the nucleons in the photon-nucleon scattering 
amplitude, which is a small correction of G(e/mjv)) where e is the average nuclear binding 
energy and is the mass of the nucleon. The effect of the off-shellness in nuclear DVCS 



was considered and estimated in Refs. 



22 



23). 



7 



The initial and final states of the active nucleon are 

\p N ) = \a(l + U)Pt,k± - <*-^,\) , 

\p' N ) = \a'(l - £a)P+ k' ± + a'^, A) . (10) 

The LC fraction and the transverse momentum of the active nucleon are found from the 
conservation of the light-cone energy-momentum in the elementary 7*iV — > ^yN vertex, 

a = - —a - — « a - 2£ A , 

A = k ± + ^—^A ± ^k ± + (l-a)A ± . (11) 
1 - 4a 

In the above equations, the approximate relations hold after one neglects £a compared to 
unity. 

The function is the overlap between the initial and final nuclear LC wave functions, 

p%(a',k' ± , X\a, k x , A) = f^/ ^^ j PnW, k'±, A)<M«, k±, A) 
1 doL;d 2 k\; 1 ' 



E / II ^J^ ^ + E Q J - *) Ifaflfe + E MI<M^ frit, A t )| 2 
<t>* N {a\ k' ± , A)0Ar(a, k±, A) . (12) 



The last line is an approximation valid for sufficiently large nuclei, when the effects associated 
with the motion of the center of mass of the nucleus (taken into account by the 5 functions) 
can be safely neglected. Note that the helicity conservation requires that the helicity of the 
active nucleon be the same in the initial and in the final state. 

The matrix element in Eq. (Q can be evaluated by making a transverse boost to the sym- 
metric frame of the active nucleon {2]]]. In that frame, one can use the standard definition, 

-i J d i xe-^ x (p' N \T{J^x)r(0)}\ PN ) = H^^ N ,t,Q 2 ), (13) 

where is the DVCS amplitude for the bound nucleon. The skewedness £jy is determined 
with the respect to the active nucleon, 

fr=7^- = 7TXiT P ' (14) 
Ap N -q (1 + f A )a - £a 



S 



where pn = (pn+Pn)/^- Therefore, we obtain the connection between the DVCS amplitudes 
for the nuclear target and for the bound nucleon, 

da d 2 k 

■ I : 

N X 

To the leading twist accuracy, the DVCS amplitude for the bound nucleon is parametrized 
in terms of four nucleon GPDs, if at, En, Hn and En'- 



„ , ./ Va a' 167^ 



If 1 r 

H N u ^N,t,Q 2 ) =— + (-9±) / dxC + (x,£ N ) H N (x,£ N ,t)u(p' N )>y + u(p N ) 
%P N J-i L 

ia+ x Ax 1 

+ E N (x,£ N ,t)u(p' N )— u(p N ) +..., (16) 

ZrriN -I 

where . . . denotes the contribution of the GPDs and En- The tensor is defined in 
the boosted frame (the symmetric frame of the active nucleon) and, to a good accuracy, is 
equal to entering Eq. ([5]), 

q-PN [q-PNr vpn-pnY 

, r -^- PV + o(^,^), (17, 

where the vectors q and Pn refer to the boosted frame; fn 2 N = m 2 N — t/4, and Ra is the 
nuclear radius. In the derivation of Eq. (fTT|) we used the fact that the transverse boost 
to the symmetric frame of the active nucleon has not changed the plus-component of the 
vectors and that the typical (transverse) momenta of nucleons in a nucleus, \pn±\ ~ ^/Ra, 
are small compared to the virtuality Q 2 . 

Using the fact that the helicities of the bound nucleon in the initial and final states are 
the same and making a natural assumption that p 1 ^ is the same for the A = ±1 helicities, 
we observe that the nucleon GPDs H and E do not contribute to Eq. (fl5|) . which is a 



consequence of the light-cone spinor algebra (see, e.g., Ref. |2l|). In addition, since we are 
interested in the kinematics, where the values of xb and £/v are small, the contribution of the 
GPDs E, which enters Eq. ( 1151) with the prefactor £ N , can be safely neglected. Therefore, 
we have that the DVCS amplitude for the bound nucleon reads (keeping in mind the equal 
helicities of the initial and final nucleon) 



Htf(Z N , t, Q 2 ) = -iff \l-e N n N (^N, t, Q 2 ) w -gT H n (£n, t, Q 2 ) , (18) 
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where 1~Ln is the CFF of the bound nucleoli. Thus, we obtain our final relation between 
the CFF of the nuclear target in the impulse approximation, Ti^ , and that of the bound 
nucleon, 

da d 2 k 



EE 

N X 



y/a a' 167T 3 



Pa (<*', k '±> M a > k ^ X) 7~Ln({;n, t, Q 2 



(19) 



It is important to point out that the integration over a (longitudinal convolution) and k± 
(transverse convolution) takes into account the effect of the motion of the bound nucleons in 
the target (Fermi motion effect). The Fermi motion effect in DVCS on nuclear targets in the 



form of longitudinal convolution was also considered in Refs. 21 



24 



25 



27|. Both the 



longitudinal and transverse convolutions along with the modifications of the bound nucleon 



22 



23] 



GPDs, which depend on k±, were considered in Refs. 

To interpret the function p A and to fix its normalization, it is useful to consider the 
electromagnetic form factor of a spin-0 nucleus, F^ m ', which is defined as the matrix element 
of the operator of the electromagnetic current, 



[p A \no)\PA) = 2P^Fr-(t) 



(20) 



Using the LC formalism just presented, we consider the plus-component of Eq. (12011 and 
obtain 

da d 2 k 



2PJFT-(*) = EE / 

AT \ J 



N A 



y/a a' 167T 3 



Pa («'> M a > k ±i A) (p' N \ J + {0)\pn) 



(21) 



In the reference frame that we work in, the momentum transfer A is predominantly trans- 
verse at small xb [see Eq. (T5])]. Therefore, the nucleon matrix element for the same nucleon 
helicities is (predominantly) proportional to the Dirac nucleon form factor, Fu\r(t), 

( P ' n \J + (0)\Pn) « u(p' N h + u(p N )F 1N (t) » 2p + N F 1N (t) = 2^P+F 1N (t) . (22) 



Therefore, 



W) = EWrE / ^^P^'A 1 X\a } k ±} X) = J2FA(t)F 1N (t), (23) 
n ^ j y a a N 

where we have introduced the nuclear form factor associated with the distribution of nucleons 
in the nucleus (associated with the nuclear density), 



F A (t) 



£a f da d 2 k± 
&v J \faa' 167r 3 

10 



p A (a', k' ± , X\a, k±, A) 



(24) 



As follows from Eq. (1231) . F A (t) is normalized to unity [-Fa(O) = 1]. This condition also fixes 
the normalization of the nuclear LC wave function, 

A A 

At small ib, the effect of the Fermi motion can be safely neglected (see, e.g., Ref. 28]), 
and, as a consequence, Eq. (1191) can be significantly simplified as follows. The function 
is peaked around a ~ 1/A. Thus, if one neglects the Fermi motion of the bound nucleon, 
one evaluates £jv at a = l/A (where, for brevity, we shall use the same notation), 

6v = frr(a = 1/A) = rTT-T^ 7 « ^ - (26) 

Therefore, neglecting the Fermi motion and using Eq. ff2M . Eq. ( fl9|) can be written in the 
following simplified and approximate form: 

£_ 

A ' AT 

As a number of nucleons, scales as A 2 , which is a natural scaling of the nuclear CFF 29]. 



The inclusion of the Fermi motion effect and the effect associated with non-nucleon degrees 
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30] 



of the freedom in the nucleus modifies this intuitive scaling 

The next important step is the conversion of the relation between nucleus and nucleon 
CFFs [Eq. ( 12 7p ] into a similar relation between the corresponding GPDs. To the leading 
twist accuracy and to the leading order in the strong coupling constant, 

H A (U,t) = [ dxH A (x,U,t)( ^— + ' 



i \x-£ A + ie x + £ A -ie 

n N (£ N ,t) = [ dx N H N (x N ,£ N ,t)( l —— + — t -V (28) 

J-i \x N -£ N + ie x N + £ N -ieJ 

The relevant quark LC fractions and momenta of the active nucleon and the target nucleus 

are presented in Fig. [31 Figure [3](a) represents the generic handbag approximation for DVCS 

on a nuclear target, which expresses the CFF 7i A in terms of the nuclear GPD H A and which 

corresponds to the first line of Eq. (T281 . 

At the same time, can be expressed in terms of the nucleon CFF Hn [see Eq. ( 1271) 

and Fig. 13(b)]. In this case, the nucleon GPD depends on the LC fractions £zv defined by 

Eq. (TbT| and on ijv, which is defined with respect to the active nucleon, 

k ' Tl X 
XN = ~ = 7— -T-, T , 29 

p N -n a(l + £ A )-£ A 
11 



7*(g) 7 7* (9) 7 




FIG. 3: The handbag mechanism for DVCS on a nuclear target, (a) The generic representation of 
nuclear GPDs. (b) A more detailed representation of the same quantity in terms of bound nucleon 
GPDs. Shown are relevant quark light-cone fractions and momenta of the active nucleon and the 
target nucleus. 



where k = (k+k')/2 and k and k! are the momenta of the initial and final lepton, respectively. 
A useful consequence of Eq. (1291) is the proportionality of the LC fractions xn and x: 

^ = f. (30) 

This relation allows us to find the LC fractions of the interacting quark in Fig. EJ^b), which 
are equal to x + = (£a/£n)(%n +6v) and x — £a = (£,a/£,n)(xn — Cat), respectively. Since 
the absolute value of xjy cannot exceed unity, we find that 

N<£-i. < 31 > 

Note that the limit \x\ < 1/A is standard for the approximation, when the nucleus consists 
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of A stationary nucleons. Using Eq. (Fl9|) and the second line of Eq. (1281) . we obtain 

x [ dx N H N (x N ,£ N ,t) ( \ — — H — -J 

x/ dzfT^xjv,^,*) + ■ ( 32 ) 

Recalling the first line of Eq. (1281) and the limits of integration over x [Eq. (I3T1) ]. we obtain 
the desired relation between the nuclear GPD in the impulse approximation, and the 

nucleon GPD: 

^(ai,Q 2 ) = VV / -^=^p^a',k' ± ,X\a,k ± ,X)H N (x N ^ N ,t,Q 2 ). (33) 

We would like to note that Eq. ([3"3"j) could also be derived starting directly from the definition 
of the nuclear GPD as the matrix element between nuclear states and applying the LC 
formalism for the nuclear states, as we did for the DVCS amplitude above. 

Equation (133]) is derived for the nuclear (nucleon) GPDs, which are sums of quark GPDs 
weighted with the quark electric charge squared. Certainly, the relation between the nuclear 
and nucleon GPDs holds for individual parton flavors (quarks and gluons): 

Hf\Ut,Q 2 ) = J2J2 I P~! jh Pa K %, k±, A) W N {x N , 6v, t, Q 2 ) , (34) 

where j is the parton flavor. 

As we have already explained, the Fermi motion effect can be safely neglected at large 
energies [see Eq. (|2"T)) ]. In this case, Eq. (I34|) can be simplified and written in the following 
form: 

Hf\x, U, t, Q 2 ) w ^ H n(*n, tN, t, Q 2 ) . (35) 

B. Double scattering correction 

The graph in Fig. [2(b) describes the contribution to DVCS on a nuclear target, when the 
interaction involves two nucleons of the target. This graph gives the leading contribution 

13 




FIG. 4: Double rescattering correction to DVCS on a nuclear target, (a) The shadowing correction 
in terms of the j*NN — > ^yNN amplitude, (b) An approximation based on the assumption that 
the shadowing correction can be expressed in terms of DVCS on a Pomeron, j*JP — > 'jJP. Also 
shown are the relevant light-cone momentum fractions. 



to nuclear shadowing. Details of the kinematics of the graph in Fig. E](b) are presented in 
Fig.H 

Using the LC formalism, we obtain the following expression for the contribution of the 
graph in Fig. [2](b): 

rr(b)^ _■ f ,4 -iq-xST f ^1 da l da 1 d 2 k ±1 d 2 k ±2 

x pf (aW2,k' ±1 ,k' ±2 \a h a 2 ,k ±1 ,k ±2 ) (p' 1 p' 2 \T{J^x)r(0)}\p 1 p 2 ) , (36) 

where J^pairs denotes the sum over the pairs of the active nucleons with momenta p\ and p 2 
in the initial state and with momenta p[ and p 2 in the final state. Each state is characterized 
by the corresponding LC fractions and transverse momenta: 

b'i, 2 ) = l«i, 2 (l - U)Pt k'±i,2 + «i,2^) • (37) 
The LC fractions and the transverse momenta of the active nucleons are related by the 
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conservation of the LC energy-momentum [see also Eq. (TTTT) ] : 

a[ + a' 2 = ai + a 2 — 2£a , 
k' ±l + k' ±2 = k ±1 + k ±2 + A ± , 



(38) 



where we have neglected the factors £4 and a^ 2 compared to unity. 

For brevity, we shall not show explicitly the nucleon helicities keeping in mind that the 
interaction does not change the helicity of the nucleons. The function is given by the 
following overlap of the nuclear LC wave functions: 



Pa («i a 2) fc±i> k' ±2 \a!, a 2 , kj_i, k ±2 ) = 4>* N (a[, k' ±1 )(p N (ai, k_Li)(j)* N (a 2 , k'_ L2 )4> N (a2, k± 2 ) 



II «; - 1) 16vr 3 5(^ k ±j ) k' A 

i=3 j=l j=l 



Li) 



0jv( a i> fc±i)0Jv(ai, k±i)(j)* N (a 2} k' ±2 )<j) N (a2, k 



L2 



(39) 



Equation (1361) is a general expression corresponding to the graph in Fig. [21(b) and to the 
graph in Fig. 0(a). To proceed with the derivation, we need to model the multiparticle 
matrix element {v'ip' 2 \T{J^{x)J v (0)}|pip 2 ) . Our model for the (p'iP2|7 1 {^ A '(^)^ iy (0)}|piP2) 
matrix element is based on the studies of hard inclusive diffraction in DIS on the proton 



31 



32 



33. 



34|, which we shall briefly review in the 



at HERA in the reaction ep — > eXp 
following. 

The diffractive DIS ep — > eXp reaction is presented in Fig. [5j The ep — *> eXp cross 




X(M X ) 



FIG. 5: Diffractive DIS on the proton. 
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section is expressed in terms of the diffractive structure functions F®^ and F®^ as 



d'< 2ira 2 em 



dx jpdtdx pdQ 2 xbQ 4 



(1 + (1 - y) 2 )F 2 D(4) (x B , Q 2 , xjp, t) - y 2 F^ A \x B , Q 2 , xjp, t) 

(40) 

where a em is the fine-structure constant and y — (p-q)/(p- k) is the fractional energy loss of 



the incoming lepton. The variables t, xjp, and (3 are characteristic for diffractive processes, 

t = (p'-p) 2 , 

q-(p- p') M\ + Q 2 



Xjp 



q-p W 2 + Q 2 ' 

x_ = q 2 g_ 2 

x F 2q-{p- p') ~ Q 2 + M 2 X 

where Mx is the invariant mass of the diffractively produced final state and W 2 = (q + p) 2 - 
The variable xjp is the fraction of the proton LC momentum lost in the diffractive scattering 
(the LC fraction carried by the Pomeron); f3 is the LC momentum carried by the interacting 
quark (parton). As follows from the definition of xjp, the minimal value of xjp is equal to 
Bjorken xb, which corresponds to Mx = 0. Typically, the contribution of Fj ^ is neglected 
because of its smallness and because of the kinematic suppression by the y 2 factor. 

One of the main physics results of HERA is the observation that hard diffraction in DIS 
constitutes a large part (10-15%) of all DIS events and that hard diffraction in DIS is a lead- 
ing twist phenomenon, that is, that the diffractive structure function F®^ approximately 
scales (i.e., it only weakly - logarithmically - depends o n Q 2 ). 



The factorization theorem for hard diffraction in DIS 35| states that, at given fixed t and 
xjp, the diffractive structure function F®^ can be written as convolution of hard scattering 
coefficient function Cj with the universal diffractive parton distributions ff^ (J is the 
parton flavor): 

F 2 ^(*,Q 2 ,W) = — £ f ^C,(^-,Q 2 )f^(P',Q 2 ,xjp,t). (42) 

It is a phenomenological observation, which follows from the QCD analysis of the HERA 
data on inclusive diffraction, that the diffractive PDFs ff^ can be written as a product 
of the Pomeron flux, fjp/ p , the parton distribution function of the Pomeron, fj/jp, and the 
factor describing the t dependence, 

ff (4) (/3, Q 2 , xjp, t) = fjp /p (xjp)f j/]P ((3, Q 2 )B diS e B ^ 1 . (43) 
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In Eq. (|43|) . we neglected the contribution of the subleading (Reggeon) exchange, which is 



not important in the considered kinematics. The Pomeron flux has the following form 



33, 



34) 



fjp/p(xp) = I dt Ajp 2ap{t) _ x , a F (t) = «jp(0) + a'jpt , (44) 



1 GeV 2 X p 



where t min » -m 2 N x 2 B « 0; fijp = 5.5 GeV~ 2 , aj>(0) = 1.111 (Fit B of Ref. |33j), and 



a 



ip 



0.06 GeV 2 . The coefficient Ap is found from the condition xjpfjp/ p (x]p) = 1 at 



X]P = 0.003. 



global fits to the HERA data on hard 
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32 



33. 



using the QCD factor- 



The PDFs of the Pomeron, fj/jp, are found from 
diffraction taken by the ZEUS and HI experiments 
ization theorem [Eq. (142j) ]. One of the main results of such fits is that the gluon diffractive 
PDF is much larger than the quark diffractive PDFs. 

The t dependence of hard inclusive diffraction at HERA was recently measured by the 
HI collaboration using the forward proton spectrometer, which allows to detect the final 



proton 34(. In the kinematics of the experiment, the data was well described by the simple 
exponential form [Eq. (1431) ] with the slope B diS ps 6 GeV -2 . (Note that has the 

dimension GeV -2 .) 

Our model for the (p'iP2\T{J^{x)J v ($)}\p\p2) matrix element is based on the observation 
that, in the considered kinematics, the interaction of the active nucleons with the virtual and 
real photons has a diffractive character and, hence, proceeds via the t-channel exchange with 
the vacuum quantum numbers (i.e. the Pomeron). The model is schematically presented in 
Fig.[4](b). The space-time picture of the process is the following. Nucleon 1 with longitudinal 
momentum fraction a\ emits a Pomeron with momentum fraction a\Xp. The virtual photon 
undergoes DVCS on that Pomeron, producing a real photon and a Pomeron with the LC 
fraction a\Xjp — 2£a, which is absorbed by nucleon 2. Note that while the skewedness £a is 
fixed by the external kinematics, the variable xjp is integrated over since it is related to the 
LC fractions of the active nucleons, 

a[ = cti — aiXjp , 

a' 2 = a 2 + aiXjp - 2^ A . (45) 

The variable xp has a clear physical interpretation: it is the fraction of the LC momentum 
of the nucleon carried by the Pomeron (see the previous discussion of diffraction in DIS). 
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Whereas xjp is the relevant variable for the Pomeron emitted by nucleoli 1, for the Pomeron 
emitted by nucleon 2, the relevant fraction is 

= ~ xp ~ 2&v • (46) 

a 2 a 2 

Based on this discussion, our model for (p'iP 2 \T{J^ l (x)J u (0)}\pip 2 ) reads 
-i J d i xe-^ x {p' 1 p' 2 \T{J^x)rm\PiP2) 

= -(2n)k v 16nB diS 4> 1P / N (x]p)(f) 1 p/ N (x ] p - 2^ N )—H^(^ 1P ,t,Q 2 ) , (47) 

Xjp 

where k v = (l — irj) 2 / (1 +i] 2 ), rj « 7r/2(«jp(0) — 1) ~ 0.17 is the ratio of the real to imaginary 



parts of the 7*iV — > XN diffractive amplitude [111. Il2l. Il3[] . (f)p/N is the probability amplitude 



of emitting a Pomeron off the nucleon, and Hv, is the DVCS amplitude on the Pomeron. 



In our analysis, we take 4>ip/n(xip) = a/ fp/ P (xip), where the Pomeron flux is defined by 
Eq. (jUj) . The DVCS amplitude on the Pomeron, Hp , is modeled by using the PDFs of the 
Pomeron, fj/p, which enter Eq. (143 p . The t dependence of Hp is given by the factor e BdiBt . 
The skewedness £jp is defined with respect to the Pomeron [compare to Eq. (111]) ]. 

i* = j?—= U p , (48) 
4 pip ■ q ot\Xp — £ A 

where pip = {pip + p'ip)/2 with pjp and p'jp the momenta of the Pomerons emitted by nucleon 
1 and nucleon 2, respectively. 

A few words are in order about the remaining factors in Eq. (l4"Tj) . The factor of 2tt comes 
from the standard definition of the connection between the Compton scattering amplitude 
and the structure functions. The factor of 16n is specific for diffraction and has its origin 

n 

in the optical theorem (see, e.g., Ref. [111]). Note also the overall minus sign, which is a 
consequence of the fact that the considered matrix element is essentially a product of two 
scattering amplitudes, which are predominantly imaginary at high-energies. 

To implement Eq. (j4"T|) in Eq. (136]) . we insert the following identity in Eq. (136]) : 

l = y da' 2 dxp 5(a' 2 — a 2 — a 2 {xp — 2£ N )) 5(a' 1 — ot\ + aiXp) ot\ . (49) 

Inserting Eq. (jlTj) in Eq. (136]) . we obtain 

crx f f da' 1 da' 2 d 2 k'_ Ll daida 2 d 2 k ± id 2 k ± 2 f ' 1 , r/ , / „^ u 

x 8(a' 1 - «i + aiXj P )ai0^(a' 1 , ^^(ai, ^±i)0at(« 2) k' ±2 )(j) N (a>2, k ±2 ) 
x k v (32n 2 )B diS (p ] p /N (x 1 p)4> I p /N (x I p - 2£ N )^-^-Hp U (£p,t,Q 2 ) , (50) 
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where = max{x7v, 2£tv}. The limits of integration over xp deserve a comment. The 
lower limit of integration is the simultaneous requirement that the Pomeron LC fraction in 
Eq. (j46p is non-negative [see also Fig. H^b)] and that the Pomeron LC fraction is larger than 
the LC fraction of the active quark, xjp > x^. The upper limit of integration is the standard 
condition on the produced diffractive masses, which can be cast in the form xp < 0.1. 

In addition, in Eq. (1501) we made an assumption that multiple interactions with the 
target nucleons lead only to the attenuation of v and do not introduce an additional 
imaginary contribution. This amounts to taking the real part of the expression describing 
the interaction with two nucleons of the target. 

For the comparison with the predictions of the LT theory of nuclear shadowing for nuclear 
PDFs and for the convenience of numerical calculations, we evaluate the overlap of the 
nuclear LC wave functions in Eq. (150]) in the coordinate space. The Fourier transform of 
the nuclear LC wave function reads 

<M<*, kx) = [ dzd 2 be imNaz+l ^<f) N (z, b) . (51) 



The normalization of the LC wave function in the momentum space [Eq. ( 1251) ] fixes the 
normalization of the wave function in the coordinate space, 

J dzd 2 b\<p N (z,b)\ 2 = J dzd 2 bp A (z,b) = 1, (52) 

where pa{z, b) is the nuclear density. We have used that a ~ 1/ A. In our numerical analysis, 
we used a two-parameter Fermi form for pa{z, b) 36]. 



Thus, substituting Eq. floTj) into Eq. (jSTfl) . using the approximation 

°' (53) 

and integrating over the light-cone fractions and the transverse momenta, we obtain our 
final expression for H^^: 

^ = _ a(a- i) |v 167r£?diff ft e j f d ^ e iA X r b r dzi r dz2 f°- dxpPA (b,z 1 )p A (b,z 2 ) 

x k v e- im - Z2 (^- 2 ^) +m - Zl ^0 JP/ ^(a; J p)0 J p /J v(a; J p - 2^ N )}—H^^ P , t min , Q 2 ) , (54) 

where we have used that J^pairs = A(A — l)/2. Note that to perform the Fourier transform, 
we neglected the weak t dependence of Hp compared to the rapid t dependence of the 
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nuclear distribution and, hence, evaluated Hp at the minimal momentum transfer t min m 
—m 2 N x 2 B ~ 0. We also introduced the z 2 > z\ ordering to reflect the space-time evolution of 
the 7*iV7V — > jNN scattering (see also, e.g., Ref. [371] ) . 

Equation (|54p can be turned into the relation between the nuclear GPD and GPD of the 
Pomeron, quite similarly to the corresponding derivation in the previous section. The DVCS 
amplitude on the Pomeron, Hp ', is expressed in terms of the CFF of the Pomeron, Hjp, as 

H^^jp, t, Q 2 ) « -gfHrtejp, t, Q 2 ) , (55) 

where we neglected the same terms as in Eq. (|17p . Therefore, for the contribution of the 
graph in Fig. [2](b) to the nuclear CFF we obtain 

n (b) = _ A (A- l) |v 167rj g d . ff ^ e | f d ^ e ^ rb f°° d Zl f°° dz 2 dxjpp A (b,z 1 )p A (b,z 2 ) 

^ kA ^ J Joo Jzi Jxf, in 

x kr, e-^^-^+^^^jp/^xjp^jp/Nixjp - 2^ n )}—Hjp^jp, t min , Q 2 ) . (56) 

) Xjp 

To the leading twist accuracy and to the leading order in the strong coupling constant, 
Tip can be expressed in terms of the GPD of the Pomeron, Hp, as 

H F (( F , t) = /' 4 W , t) (— ^ + ) . (57) 

Using the same argument that led to Eq. (1501 . we find that 

X X X]\f / r o \ 

T ~ T ~ ~T ' ^ > 

C,P 

where x parametrizes the interacting quark LC fractions in the graph in Fig. El^a). Those 
fractions are equal toi + ^ = (£4/ '£jp)(x' + and x — £a= (£a/6p)( x ' — £zp), respectively. 
Since \x'\ < 1, we find that \x\ < £a/£p- Thus, substituting Eq. (156"j) into the first line 
of Eq. (jUJ), changing the integration variable from 2 to x' according to Eq. fl58l) . recalling 
Eq. (1571) . and noticing that the ensuing relation holds not only for the DVCS amplitude 
written to the leading order in the strong coupling constant, but also for individual parton 
flavors, we obtain the contribution of the graph in Fig. [2(b) to the nuclear GPD of flavor j, 



H 3 />, as 



3(b) 

Hf\x,U,t,Q 2 ) = - Ai - A ~ l) ^lQixB A ^e{ [ d 2 be iK ^ b f° dz x f°° dz 2 f' dx F 

X PA (b, Zx)p A (b, Z 2 ) k v e -im N z 2 (x 1P -2^ N )+im N z 1 x F 

x 4>p/n{xip)4>ip/n{x]p - 26v) } — H j jp{-f-x N , £jp, t min , Q 2 ) . (59) 

J Xjp t, N 
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The GPD of the Pomeron, Hjp, is modeled by using the PDFs of the Pomeron, fj/jp- In 
our numerical analysis, we used the model of GPDs in which it is assumed that the effect of 
skewedness in GPDs can be neglected at the initial evolution scale. This model corresponds 
to the double distribution model 38] with a 5-function-like profile [39]. The details are given 
in Sec. HV] 



C. Quasi-eikonal approximation for multiple rescatterings and the final expression 
for nuclear PDFs 

To evaluate the contribution of the graph in Fig. Wis), we use the following high-energy 
(small xb) space-time development of the process. The virtual photon diffractively interacts 
with nucleon 1 and produces a certain diffractive state X characterized by xjp (diffractive 
mass Mx)- The produced state is then assumed to elastically scatter on A — 2 nucleons of 
the target. Finally, the last interaction of the state X with nucleon 2 produces the final real 
photon. This picture of multiple rescattering at high-energy corresponds to the quasi-eikonal 
approximation for the graph in Fig. [2](c) and higher rescattering terms. The quasi-eikonal 
approximation was used in the evaluation of nuclear PDFs in the framework of the leading 

nnn 

twist theory of nuclear shadowing pjj, [12|, H3J and in the evaluation of the DVCS amplitude 



on nuclei in the framework of Generalized vector meson dominance model [21 

Within the quasi-eikonal approximation, the multiple interactions can be summed and 
can be cast in the form of the eikonal attenuation factor, T, 

T = e -^a-ivK a (xB,Q 2 ) f z J dz' PA (b,z') ^ ^ 

where al s is the effective cross section, which determines the strength of the rescattering of 
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off 

the state X off the nucleons. This cross section is defined as 

v 3 efi {x B ,Q 2 ) = . 2 / 67r f dlff , n2 , [ dxjpPfjp/pix^fj/jpi^Q 2 ), (61) 

[i+r] )XBjj/N{XB, W ) Jx B 

where fj/N is the usual parton PDF of the nucleon. For a given flavor j, a J eS is proportional 
to the probability of diffraction relative to the total probability of the interaction. As an 
example, we present cr J c{i as a function of x B at fixed Q 2 = 2.5 GeV 2 for the u quark and 
gluon flavors in Fig. [61 
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FIG. 6: The effective cross section a J eS [see Eq. ([51])] for the u quarks and gluons as a function of 
Bjorken xb and at fixed Q 2 = 2.5 GeV 2 . 



Thus, collecting all contributions to the nuclear GPD H J A , 



Hi(x, U t, Q 2 ) = Hf> + H j } h) + HT + ..., 
we obtain our final expression for flavor j GPD of a heavy spinless nucleus: 
Hfc, U, W 2 ) = t-FaV) J" H n(^n, 6v, t, Q 2 ) 



T 3(c) 



(62) 



^4-i)e 



N 



6 



j2r AA I b 



poo /"0.1 



x p A (6, «i)p A (6, z 2 ) k v e- imNZ ^ XF ~ 2 ^ +imNZlXF e-^ {1 ~ ir ^^^ 

(63) 



x (f)]p/N(x]p)4> r / N (x]p - 2£ N )\—Hjp(-^-x N ,£]p,t miri ,Q 2 



For practical applications and for a comparison to the case of a free nucleon, it is convenient 
to simultaneously rescale the LC fraction x and the nuclear GPDs in the left-hand side of 
Eq. (El: 



H J A (x,Ut } Q 2 



N 



6 



H A (x N ,U,t,Q 2 ) , 



(64) 



(where the rescaling of the nuclear GPD is necessary to preserve sum rules involving the 
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nuclear GPD). Then, our master equation for the nuclear GPD becomes 
H A (x N , U, W 2 ) = F A (t) H j N {x N , £ N , t, Q 2 ) 

N 

A(A-l) - r - - r °° r °° r0A 



16nB dlS ^e\ / d 2 be iK ^ / dz x \ dz 2 / dx P 

^ J Joo Jzi Jx™ in 



2 

X PA (b, Zx) PA (b, Z 2 ) k v e -^MZ2(x P -2 M+ irn N z lXFe -^l-ivH B ^B 1 Q 2 )f z J dz> PA {b,z') 

x 4>iP/N(xip)<f>ip/N(xip - 2£at) \ — H J jp(-^-x N ^ip,t min ,Q 2 ) . (65) 

} X]P t,N 

As we explained above, we neglected the Fermi motion effect in the first term in Eq. fl65l) . 
If necessary, the Fermi motion effect can be restored by replacing the first term in Eq. (1651) 
by the right-hand side of Eq 



III. NUCLEAR GPDS IN THE U -> LIMIT AND THE SPACIAL IMAGE OF 
NUCLEAR SHADOWING 

In the forward limit, nuclear GPDs reduce to nuclear PDFs, 

H A (x,0,0,Q 2 )=f, /A (x,Q 2 ). (66) 
Taking the £ A = t = limit in Eq. fl65l) . we obtain 
H A (x B ,0,0,Q 2 ) = J2fiMx B ,Q 2 ) 

N 

A(A — 1) r f - [°° f°° f 0A 

167r£> diff 3?e< / d 2 b / dz\ \ dz 2 / dx P 



OO J Z\ J Xq 

A , 



2 

x p A (b, z 1 )p A (b, z 2 ) k v e m ^(zi-z2) e -^-^K H (xB,Q*)f z J dz' PA{ b,z') 
x fjP/ P (x P )—f j/P {P = — ,Q 2 )) • (67) 

Xjp Xjp J 

Here we used the fact that, in the £ A — > limit, £at, £jp, £ m in — ► and $,ipXn/C,n = x' — > 
(3 = Xb/x p . The obtained expression for the nuclear PDF fj/ A as a forward limit of the 
nuclear GPD coincides with the direct calculation of fj/ A in the framework of the leading 



twist theory of nuclear shadowing [ill, [12J, ll3j; that is, our master equation [Eq. fl65l) ] has 
the correct (consistent) forward limit. 

Next let us consider the £ A — > limit (i.e., the limit when the momentum transfer t is 
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purely transverse, t = —A 2 ± ). Taking the £4 — > limit in Eq. (|65|) . we obtain 
H A (x N , 0, t, Q 2 ) = F A (t) H 3 N {x N , 0, t, Q 2 ) 

N 

~ A ^ 9 ~ ~ 167rg diff ^e{ / d 2 be iK ^ j dzi [ dz 2 f dx P 

x pa(&, *i W(&, 22) k n e^^-^e-^ 1 '^^'^ % dz ' PA ^ 
xMpM4p(^ 2 )}. (68) 

Xjp Xjp ) 

Again, we used the fact that, in the £4 — > limit, &v> 6p> *min — > 0, £jp£tv/£tv xn/xjp, 
and Hpi^XNi i min , Q 2 ) — >■ fj/p(xN /xjp,Q 2 ). Note also the lower limit of integration 
over xjp, x}p m = xjy. Since the t dependence of the nuclear form factor, F A (t), is much faster 
than that of the nucleon GPD H 3 n {xni 0, t, Q 2 ), the latter can be evaluated at t — (i.e., in 
the forward limit). Then, Eq. ( 1681) becomes 



H J A (x N , 0, t, Q 2 ) = F A (t) fj/N(x N , Q 2 



N 

A ^ A ~ ^ 167r5 difT Ke| J d 2 be iK ^ J dz x [ dz 2 [ dxjp 

A, 



OO J Z\ J Xpj 



2 

x p A (b, Zl )p A (b, z 2 ) k v e ^ F (^2) e ~4(i-ivH^B,Q 2 ) f% dz> PA $,*') 

xfjp/ P M^f j/P (^^ 2 )}- (69) 

X]p Xjp ) 

In the case of nucleon GPDs, the interpretation of GPDs in the £ — > limit is given in 
the impact parameter representation, where the GPDs have the meaning of the probability 
densities [l]. We shall also analyse our nuclear GPDs in the £ A — > limit in the impact 
parameter space. To this end, we introduce the nuclear GPD in the impact parameter space, 

H A (x,0,b,Q 2 ) = J |^e-^^(x,0,t = -Ai,g 2 ). (70) 



The Fourier transform of Eq. ( 1691) gives 

H A (x N , 0, b, Q 2 ) = T A (b) ^ fj/N(x N , Q 

N 

A(A - 1 



N 

Uroo poo rO.l 

dz\ j dz 2 / dxjp 
OO J Z\ J Xpf 



2 

X PA (b, Zl ) PA {b, Z 2 ) j™»xr{zi-zz) e -^-iri)<r>^ B ,Q*)S% dz' PA (b,z>) 

xfjp/pM— /,/ip(— ,Q 2 )), (71) 

Xip Xjp J 
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where Ta(6) = f dzpA^b, z) and pA(b, z) is the nuclear density [see Eq. (1521) ]. It is important 
to note that the nuclear GPD H J A (xN,0,b,Q 2 ) given by Eq. flTTj) is nothing else but the 
impact-parameter-dependent nuclear PDF introduced and discussed in the framework of 
the leading twist nuclear shadowing 

In Eq. ( 1711) . the first term is the Born approximation to H J A corresponding to the graph in 
Fig. [2](a); the second term is the nuclear shadowing correction corresponding to the graphs 
in Figs. Mjo) and [2](c) and to higher rescattering terms not shown in Fig. [2l We quantify 
the magnitude of the nuclear shadowing correction by considering the ratio 

TA{b)2^ N fj/N{XN,Q 2 ) 

where the numerator is given by Eq. (ITT]) . In the absence of nuclear shadowing, 
R? (xN,b,Q 2 ) = 1. The ratio R^xjy, b, Q 2 ) for the nucleus of 208 Pb as a function of xjy 
and b at fixed Q 2 = 2.5 GeV 2 is presented in Fig. [JJ In the figure, the top panel corresponds 
to u quarks; the bottom panel corresponds to gluons. 

Essentially, Fig. [7J presents the impact parameter dependence of nuclear shadowing, or 
the spacial image of nuclear shadowing. Several features of Fig. [JJ deserve mentioning. First, 
the amount of nuclear shadowing [the deviation of Ri(xN, b, Q 2 ) from unity] increases as one 
decreases x^ and b. Second, nuclear shadowing for gluons is larger than for quarks. For 
instance, at xn = 10 -5 and at b = 0, R 9 = 0.073, but R g = 0.23. Third, nuclear shadowing 
induces non-trivial correlations between xn and b in the nuclear GPD H A (x,0,b,Q 2 ), even 
if such correlations are absent in the free nucleon GPD. [In Eq. (J7TJ) we neglected the x^-b 
correlations in the nucleon GPDs by neglecting the t dependence of H j n (xn, 0, t, Q 2 ).] In 
this respect, the spacial image of nuclear GPDs at small x^ is very different from the case 
of the free nucleon: Whereas the free nucleon GPDs become independent of b in the x^ —> 
limit [7], the suppression of nuclear GPDs by nuclear shadowing is strongly correlated with 
the impact parameter b. 



IV. NUCLEAR SHADOWING AND PREDICTIONS FOR NUCLEAR DVCS OB- 
SERVABLES 

It is convenient to quantify the amount of nuclear shadowing in our master expression 
for the nuclear GPD of a heavy nucleus [Eq. (|65|) ] in terms of the R^xjy^NytyQ 2 ) ratio, 
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FIG. 7: Impact parameter dependence of nuclear shadowing for 08 Pb. The graphs show the ratio 
R?(xn, b, Q 2 ) of Eq. (|72[) as a function of the LC fraction and the impact parameter b at fixed 
Q 2 = 2.5 GeV 2 . The top panel corresponds to u quarks; the bottom panel corresponds to gluons. 
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which we define as 

R j (x N £ N ,t,Q 2 ) = 



2 . H A (x N ,£ A ,t,Q 2 ) 



F A (t) J2 N H' N (x N ^ N ,t,Q 2 ) 

•OO POD pO.1 



1- n 167rff diff 3fe<! I d 2 be iA± b I dz x I dz 2 I dx 



x p A (b, z 1 )p A (b, z 2 ) k n e-^^(x^-2 ?iV ) + i miV , 1 ^ e -4(i-^)^ ff (- B ,Q 2 )/; i 2 dz' PA (b,z') 
x 4>ip/n{xip)<P]p/n{x]p - 26v)| —H 3 jp{-^-x N , 6p,*min> Q 2 ) 

) X]p t,N 

F A (t)J2 H ^N^N,t,Q 2 )) . (73) 

N J 

In the absence of nuclear shadowing (and the Fermi motion effect), Ri (xiy,^iy,t,Q 2 ) = 
1. The ratio R^(x^,^N,t,Q 2 ) is a generalization and a Fourier transform of the ratio 
Ri(x N ,b,Q 2 ) of Eq. ([72D. 

At high energies, scattering amplitudes are predominantly imaginary. As follows from 
Eq. fl28|) . to the leading twist accuracy and to the leading order in a s , the imaginary part 
of the DVCS amplitude (the CFF) reads 

SmH A (U, t, Q 2 ) = -kH a (U, U, *, Q 2 ) , (74) 

where 

h a {u, a, t) = J2 e l H A(^, a, t) . (75) 

Therefore, in our numerical analysis that follows, we shall present our predictions for 

R'^N^N,t,Q 2 ). 

In our numerical analysis, we use the model of GPDs of the free nucleon and the Pomeron, 
in which it is assumed that the effect of skewedness in GPDs can be neglected at the initial 
QCD evolution scale (Ql = 2.5 GeV 2 in our case). Then, in the xn = £/v case of interest, 
one has 

H J N (£,N,£,N,Qo) = fj/N(£,N,Qo) 

H J A^ 6p, Ql) = fj/jpfo, Ql) = Ql) • (76) 

Xp 



This model corresponds to the double distribution parameterization of GPDs 



with a 5- 



function-like profile [39[]; we shall refer to this model of the GPDs as the forward-like model. 
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FIG. 8: Nuclear shadowing for the DVCS amplitude for 208 Pb at Aj_ = 0. The plots show the 
ratio Ri(£,N,£,N,t, Q 2 ) of Eq. (|73|) as a function of xb at fixed Qq = 2.5 GeV 2 for the forward-like 
model of GPDs (solid curves). For comparison, the ratio of the usual nuclear to nucleon PDFs, 
R^(xb) = fj/A( x B, Qo)/[Afjm(xB, Qo)], is given by the dotted curves. The left panel corresponds 
to u quarks; the right panel corresponds to gluons. 



Note that the suggestion that the GPDs at small xb and at the low input scale Qo can be 



40]. 



well approximated by the usual forward PDFs was first proposed in Ref. 

It is very important to point out that the recent analysis of the high-energy HERA data 
on DVCS on the proton unambiguously indicated that the description of the data at the 
leading order accuracy requires almost no skewedness effect in the input GPDs 41]. This 
clearly favors the forward-like model of the PDFs over other small-a^ parameterizations 
(see, e.g., Ref. [42]). 

Let us first examine the R^(C,n, 6v, t, Q 2 ) ratio of Eq. (!73|) in the situation when the 
momentum transfer t is purely longitudinal, Aj_ = and t = t min w — i^mjy. Figure [8] 
presents R^^n, £jv, t, Q 2 ) for 208 Pb as a function of Bjorken xb at fixed Ql = 2.5 GeV 2 (solid 
curves). Also, for comparison with nuclear shadowing in usual nuclear PDFs, we present 
the ratio R^{xb) = fj/A{ x B,Qo)/[Afj/N{xB,Qo)] by the dotted curves [ijj]. The left panel 
corresponds to u quarks; the right panel corresponds to gluons. 

As one can see from Fig. [SX the suppression of R-'(^N,^N,t,Qo) by nuclear shadowing 
is very large and it is larger than the suppression of R^(xb,Qq) in the forward case. This 
is one of new results of this work and it comes from our model for graph in Fig. H(b). In 
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particular, we assumed that the matrix element 

(p[p' 2 \T{J^x) r(0)}\ Pl p 2 ) cc <j> P/N {x]p)<f> P , N {x P - 2£ N ) , (77) 

which leads to the dynamical enhancement of nuclear shadowing because 4> p /n{ x p~^n) ^ 
4>ip/n(xip) for xjp close to 2£at- 

We stress that our results presented in Fig. [8] have an exploratory character and 
are subject of significant theoretical uncertainties, which include our modeling of the 
(p'iP2\T{J lJ- (x)J u (0)}\piP2) matrix element, the choice of the model for the nucleon and 
Pomeron GPDs, and the extrapolation of the fits for diffractive PDFs fp to unmeasured 
kinematic regions. 

We also mention that the rapid approach of R^(^N,^N,t,Qo) to unity as %b — > 0.1 is 
driven both by the decrease of the nuclear shadowing term and by the decrease of the Born 
term driven by the nuclear form factor at t — t min ~ —x 2 B m 2 N , FA(t m - m ). 

Next we examine the ratio W(^N,^,N,t,Q 2 ) at fixed t as a function of xb- In this case, 
the transverse momentum transfer is no longer vanishing: |A^| 2 « — A^m 2 N — t. Our results 
are presented in Fig. [91 The left panel corresponds to u quarks; the right panel correspond 
to gluons. The solid curves correspond to t = —0.005 GeV 2 ; the dotted curves correspond 
to t = —0.01 GeV 2 . For comparison, the ratio Ri(£,N,£,N,tmin,Q 2 ) at t = t min is given by 
the dot-dashed curves (the same curves as in Fig. [8]). 

As one can see from Fig. O the effect of nuclear shadowing [the deviation of 
W(^N,^N,t,Q 2 ) from unity at small xb] increases with increasing \t\. This is a natural 
consequence of the fact the Born term, whose t dependence is given by F^it), decreases with 
increasing \t\ faster than the shadowing correction term. 

Next we turn to a discussion of observables measured in DVCS. In lepton-nucleus scat- 
tering, it is convenient and natural to use the invariant energy per nucleon. For our results 
presented in the following, this means that we replace £a £,n and assume that the in- 
variant energy, a/s, is given per nucleon. Results of high-energy DVCS measurements are 
usually presented in terms of the DVCS cross section at the photon level, 

\A DYC s^N,t,Q 2 )\ 2 , (78) 



d0"DVCS ™ e 2 m X B| i it- , r\l\\1 



dt Q 4 

where a em is the fine-structure constant. For the DVCS amplitude at high energies, we use 
the leading twist and leading order in a s expression [see Eqs. and (175]) ]. 

|-4 D vcs(6v,t,Q 2 )| 2 « \H A (£ N ,t,Q 2 )\ 2 « n 2 (H A (^ N ^ N ,t,Q 2 )) 2 , (79) 
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FIG. 9: Nuclear shadowing for the DVCS amplitude for 208 Pb at fixed t. The plots show the 
ratio R? (£,N,£,N,t,Q 2 ) of Eq. (|73|) as a function of xb at fixed Qq = 2.5 GeV 2 . The left panel 
corresponds to u quarks; the right panel corresponds to gluons. The solid curves correspond to 
t = —0.005 GeV 2 ; the dotted curves correspond to t = —0.01 GeV 2 . For comparison, the ratio 
RH€n, £,n, *min> Q 2 ) at t = t m \ a is given by the dot-dashed curves. 



where Ha(£n, 6v, t) = J2 q q e lH A (^,N, £,n, t) and H A (£ N , £ N , t) are given by our master equa- 
tion [Eq. fl65|) ]. Since gluons enter the DVCS amplitude at the one-loop level, we do not use 
our results for the gluon nuclear GPD in our calculations presented in the following. Note 
also that since we do our calculations at fixed Qq = 2.5 GeV 2 , we use four quark flavors. 

The DVCS process competes with the purely electromagnetic Bethe-Heitler (BH) process. 
The BH cross section at the photon level can be written in the following form (see, e.g., 



Ref. 



39|) 



do 



BH 



Treat 



2tt 



|-4bh(6v, t, Q 



2\|2 



(80) 



dt AQH(1 -y + yi/2) J 2vr VMW) 
where y is the fractional energy loss of the incoming lepton, <p is the angle between the 
lepton and hadron scattering planes, V\ and V2 are proportional to the lepton propagators, 
and |-4bh(6v, t, Q 2 )\ 2 is the BH amplitude squared, which can be expressed in terms of its 
Fourier harmonics c^ H 39| as 



-4bh(6v, t, Q 2 , <p)\ 2 = C + E c n U cos (W>) • ( 81 ) 

or the case of a spinless 



n=l 



The Fourier harmonics for a spinless target are given in Ref. 43] . 
nucleus, \Abk(£n, t, Q 2 )\ 2 oc [Fyi(t)] 2 ; see further details in Ref. 



44|. 
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FIG. 10: Nuclear DVCS and BH cross sections for 208 Pb as a function of \t\ at fixed Q 2 = 2.5 GeV 2 
and xb = 0.001. For comparison, the DVCS cross section on the proton is given by the dot-dashed 
curves. For the evaluation of the BH cross section, we used y = 0.31 (see the text). 



Figure [10] presents our predictions for da-ovcs/dt and dasn/df for 208 Pb as a function of 
\t\ at fixed Qq = 2.5 GeV 2 and xb = 0.001. In addition to the input discussed above, for 
the evaluation of the BH cross section, we used y = 0.31, which corresponds to the highest 
among discussed energy options of the future Electron-Ion Collider (EIC), -Ei ep ton = 20 GeV 
and -^nucleus = 100 GeV/nucleon Q]- Also, for comparison, we give the DVCS cross 
section on the proton in the same kinematics (dot-dashed curves). 

Several features of Fig. [TO] deserve a discussion. First, the t dependence of the BH and 
DVCS cross sections repeats the pattern of [i^(i)] 2 with several distinct minima. In the case 
of the DVCS cross section, the minima are shifted because of the presence of the shadowing 
correction. Second, at small \t\, the BH cross section is much larger than the DVCS cross 
section owing to the enhancement by the 1/t kinematics factor [see Eq. (180j) ]. As one 
increases \t\ > \t m i n \, the two cross sections become compatible. Moreover, near minima of 
the nuclear form factor, the BH cross section becomes very small and, hence, the process 
is dominated by the DVCS cross section. Therefore, the measurement of the eA — > e^A 
differential cross section at the momentum transfer t near the minima of the nuclear form 
factor will provide a clean probe of nuclear shadowing in nuclear GPDs and nuclear DVCS 
owing to the suppressed BH background and the suppressed unshadowed Born contribution 
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FIG. 11: The i-integrated nuclear DVCS and BH cross sections for 208 Pb as a function of xb at 
fixed Qq = 2.5 GeV 2 . For comparison, the DVCS cross section on the proton is given by the dot- 
dashed curves. For the evaluation of the BH cross section, we used two energy settings: yfs = 32 
GeV (the upper dashed curve) and y/s = 90 GeV (the lower dashed curve) (see the text). 



to the DVCS amplitude. 

Next we study the t-mtegrated DVCS and BH cross sections at the photon level, 

T min a* aD vcs 

C DVCS — / "t — "77 — j 

J-l GeV 2 at 

o-bh = r n dt ^ . (82) 
J-l GeV 2 dt 

Figure [TT] presents the t-integrated DVCS and BH cross sections for 208 Pb as a function 
of xb at fixed Ql = 2.5 GeV 2 . For comparison, the dot-dashed curve shows the DVCS 
cross section on the proton in the same kinematics. For the BH cross section, we give two 
curves, which correspond to two different values of the cm. lepton-nucleus energy y/s: The 
upper curve corresponds to the low-energy option for the future EIC, -Eiepton = 5 GeV and 
-E'nucieus = 50 GeV/nucleon (y/s = 32 GeV); the lower curve corresponds to the high-energy 



option with £i op ton = 20 GeV and Nucleus = 100 GeV/nucleon (y/s = 90 GeV) (18|, ll9|. 

As one see from Fig. [HI in the discussed kinematics, the BH cross section is much larger 
than the DVCS cross section for xb < 0.01 for both considered high-energy options (lower 
BH curve) and for xb < 0.05 for the low-energy option (upper BH curve). Therefore, as far 
as the t-integrated eA — ► ejA cross section is concerned, it appears rather challenging to 
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extract a small DVCS signal on the background of the dominant BH contribution. However, 
the high luminosity of the future EIC should allow one to measure the t dependence of 
the DVCS and BH cross sections, which will tremendously increase the potential to probe 
nuclear GPDs in the domain of nuclear shadowing (small xb) (see Fig. [TD]and the previous 
discussion) . 

Another possibility to study nuclear GPDs in the small xb region is given by the measure- 
ment of DVCS cross section asymmetries (with polarized lepton beams or with lepton beams 
with the opposite electric charges), which are proportional to the interference between the 
DVCS and BH amplitudes. As an example, we consider the DVCS beam-spin asymmetry, 
Am, measured with the polarized lepton beam and an unpolarized target (which is always 



the case for spin-0 nuclei that we consider' 
for Am for a spinless nuclear target reads 
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To the leading twist accuracy, the expression 
Ac 



, 8K(2 -y)ZF A (t)^mH A ^N,t,Q 2 ) sin 

AluW) = 



^bh(6v, t, Q\ <P)\ 2 + '*'*y(*> 4(l - y + y 2 /2)\^mH A ^ N , f, Q 2 )\ 2 ' 

(83) 



where K oc V t m m — i is the kinematic factor 39], Z is the nuclear charge, Qm7i A is 
the imaginary part of the nuclear DVCS amplitude given by Eqs. ( 1731) . (1751) and (|65l) . 
|*<4bh(£a, t, Q 2 , 4>)\ 2 is the square of the BH amplitude ( IHTl) . and the minus in front cor- 
responds to the electron beam. To consistently work to the leading twist accuracy, one 
should use only the leading twist contributions to Vi(<f>), Vi{§) and |^4bh| 2 in Eq. (1831) . 
However, in the kinematics that we consider, t < 0.2 GeV 2 , Q 2 = 2.5 GeV 2 and <fi = 90°, 
the higher twist corrections are either absent (the terms being proportional to cos<p) or 
numerically insignificant, so that we simply use the standard expressions for Vi((f>), Viify 
and |^bh| 2 E). 

Figure [121 presents our predictions for A LU ((p) as a function of t as fixed xb = 0.001, 
Ql = 2.5 GeV 2 , and the angle = 90°. For a comparison, the dotted curve presents A LU for 
the proton in the same kinematics. Both curves correspond to the incoming lepton fractional 
energy loss y = 0.31, which in turn corresponds to the high-energy option of the future EIC 
with -Eiepton = 20 GeV and -E nU cieus = 100 GeV/nucleon. 

Our predictions for Am for 208 Pb are rather remarkable. The sole reason for the oscilla- 
tions of Ajjj for 208 Pb is nuclear shadowing! The trend of the oscillations can be understood 
as follows. At t — t m i n , Am = because of the kinematic factor K = (resulting from 
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FIG. 12: The DVCS beam-spin asymmetry, Am{4> = 90°), for 208 Pb as a function of t at fixed 
xb = 0.001 and Qq = 2.5 GeV 2 (solid curve). For a comparison, the dotted curve presents Am 
for the proton in the same kinematics. The calculations correspond to y = 0.31. 

the vanishing |Aj_| = 0). As one slightly increases |t| > |t m m|, the kinematic factors rapidly 
increase Am (which is clearly seen for the proton), but, at the same time, the nuclear shad- 
owing correction decreases the imaginary part of the nuclear DVCS amplitude, ^rriHA- As 
a result, Am increases, but not as rapidly as for the free proton case. At some rather small 
values of t, \t\ ~ 0.01 GeV 2 (a value that can be read off the left panel of Fig. [S]), QtuTIa 
changes sign and Am goes through zero. Note that at this values of t, the nuclear form 
factor, Fji(t), is still positive. As one increases \t\ further, \^itiHa\ increases, which increases 
\Am\ (with both QmTi,A and Am being negative at this point). As |t| is increased even 
further, the nuclear form factor changes sign and makes Am positive. The asymmetry 
stays positive until ^suiKa changes sign and becomes positive again [the form factor -Fa(£) 
still being negative]. As \t\ is increased, the mechanism of the oscillations just described 
repeats itself. 

We emphasize that the oscillations of Am are caused by nuclear shadowing that has a 
weaker t dependence than that of the Born contribution [see Eq. fl65l) ]. If the shadowing 
correction in Eq. (|65p is neglected, then the t dependence of the DVCS and BH contributions 
is the same and is given by the nuclear form factor Fj^if). Then, in the expression for the 
beam spin asymmetry, A LU , the t dependence from F^t) cancels and A LU for a heavy nuclear 
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target has the same t dependence as Alu for the free proton (i.e., without the oscillations). 



V. SUMMARY AND DISCUSSION 

We generalized the leading twist theory of nuclear shadowing for usual nuclear parton 
distributions to nuclear generalized parton distributions for quarks and gluons. We estimated 
quark and gluon GPDs of spinless nuclei and found very large nuclear shadowing. 

In the limit that the momentum transfer is purely transverse, £4 = £jv — 0, after Fourier 
transform, our nuclear GPDs become impact-parameter-dependent nuclear PDFs. Nuclear 
shadowing induces non-trivial correlations between the impact parameter b and the light- 
cone fraction x. 

Using our expressions for nuclear GPDs, we made predictions for the cross section of 
deeply virtual Compton scattering on the heavy nucleus of 208 Pb at high energies (in the 
kinematics of the future EIC). We also calculated the cross section of the purely electromag- 
netic Bethe-Heitler process and addressed the issue of the extraction of the DVCS signal, 
and, hence, the extraction of information on nuclear GPDs and nuclear shadowing, from the 
measurement of the eA — > ejA process. Based on our studies, we can propose two strategies. 
First, the eA — > e^A differential cross section at the momentum transfer t near the minima 
of the nuclear form factor is dominated by the DVCS cross section, which should allow for 
a clear extraction of the latter. Second, nuclear shadowing leads to dramatic oscillations of 
the DVCS beam-spin asymmetry, A^u, as a function of t. The position of the points where 
A LU changes sign is directly related to the magnitude of nuclear shadowing. 

It is important to note that the t variations of the DVCS and BH differential cross sections 
and the DVCS beam-spin asymmetry, Am, are very rapid, with the typical frequency of 
the order of 1/R\- This certainly poses a challenge for the future experiment since a rather 
high resolution in t will be required. 

One should also note that nuclear GPDs at small x will be accessed in ultraperipheral 
nucleus- nucleus collisions at the LHC [45]. In these collisions, the involved nuclei serve as 
sources of real photons, which enables one to study photon-nuclear interactions at energies 
up to ten times larger than those achieved at HERA. Nuclear GPDs will be accessed in 



exclusive photoproduction of heavy vector mesons 



461 ] and lepton pairs [47|. 
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